Interacting classical dimers on the square lattice 
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We study a model of close-packed dimers on the square lattice with a nearest neighbor interaction 
between parallel dimers. This model corresponds to the classical limit of quantum dimer models 
[D.S. Rokhsar and S.A. Kivelson, Phys. Rev. Lett. 61, 2376 (1988)]. By means of Monte Carlo and 
Transfer Matrix calculations, we show that this system undergoes a Kosterlitz-Thouless transition 
separating a low temperature ordered phase where dimers are aligned in columns from a high 
temperature critical phase with continuously varying exponents. This is understood by constructing 
the corresponding Coulomb gas, whose coupling constant is computed numerically. We also discuss 
doped models and implications on the finite-temperature phase diagram of quantum dimer models. 



The model of lattice coverings by hard-core dimers has 
a long history in classical statistical physics, culmi- 
nating in its solution for arbitrary planar graphs 

3, the 

calculation of correlation functions j3] , and its connection 
to Ising j3| or height models With the introduction 

of quantum dimer models (QDM) by Rokhsar and Kivel- 
son 0, and later in connection with resonating valence 
bond (RVB) and fractionalization physics in 2D Q and 
3D 0, dimer models have regained interest. 

As a step towards the understanding of the finite- 
temperature phase diagram of the QDM, we here study 
its limit when kinetic terms vanish, which is a classi- 
cal dimer model with an interaction that tends to align 
neighboring dimers. The dimer-dimer correlations which 
are critical in the noninteracting model Q are found to 
to remain critical, but with continuously varying expo- 
nents, down to a finite temperature where the dimers or- 
der through a Kosterlitz-Thouless (KT) transition [loj . 
Our results are naturally understood in the Coulomb gas 
approach |0 to the associated height model H, • In- 
deed, this interacting dimer model provides one of the 
simplest realizations (together with the six-vertex and 
the XY model) of a Coulomb gas (CG) with continu- 
ously varying exponents. We further study the effect of 
doping the model with monomers, and comment on the 
consequences for the quantum case. 

Model — We study a model of interacting dimers on 
the square lattice, first introduced in the liquid crystal 
context [l^ . defined by the partition function 



Z = ^ exp 



T 



(1) 



where T is the temperature, the sum runs over all dimer 
coverings {c}, and -/V'^(z)+iV^(l l) is the number of plaque- 
ttes with parallel dimers in configuration c. We mainly 
consider the case v = —1, with aligning interactions be- 
tween dimers. Note that the model Eq. corresponds 
to the diagonal part of the QDM Hamiltonian Q ■ 

At r = oo, Eq. describes the classical dimer cov- 




FIG. 1: The four possible columnar ground states. To each 
plaquette with a pair of parallel dimers we assign an energy 
of +v (shaded plaquettes in the figure). 



ering problem, and thus displays critical correlations j^. 
At T = 0, the dimers order in columns and the ground- 
state is four-fold degenerate, as shown in Fig. ^ This 
columnar phase breaks translational and rotational sym- 
metry. 

We study Eq. (HJ with Transfer Matrix (TM) tech- 
niques and with a Monte Carlo (MC) directed loop algo- 
rithm ^3 . The MC simulations are made on N ~ L x L 
lattices with periodic boundary conditions (PBC). 

Transition to the columnar phase — Previous simula- 
tions of the QDM on the square lattice found the exis- 
tence of a plaquette phase (which breaks translational but 
not ^-rotational symmetry) Qj. In the MC simulations 
of the classical model Eq. iQJ, we also find important 
plaquette correlations on finite samples, which neverthe- 
less vanish in the thermodynamic limit [15,] . However, 
they affect the finite-size scaling analysis of the order pa- 
rameter used in Ref. For this reason we resort to 
quantities insensitive to plaquette correlations to detect 
the entrance into the columnar phase. Following Ref.Q, 
we define the dimer rotational symmetry breaking (DSB) 
and pair rotational symmetry breaking (PSB) order pa- 
rameters as: 



DSB = TV^i lAr"(-) -iV"(l)l 
PSB = N-^\N%Z.) ~ N^ilijl, 



(2) 



with A^^(— ) (A^^(l)) the number of horizontal (vertical) 
dimers in the configuration c. Both quantities show long- 
range order in a columnar phase, but vanish in a plaque- 
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FIG. 2; (Color online) Dimer symmetry breaking cumulant 
(a), Pair symmetry breaking cumulant (b), and Specific heat 
per site (c) versus temperature for different system sizes. The 
dashed lines indicate the estimate Tc = 0.65(1). 

tte phase. To locate the transition with high precision, 
we study Binder cuniulants 17], defined as Sdsb = 1 ~ 
(DSB4)/(3(DSB2)2) andSpsB = 1-(PSB4)/(3(PSB')2); 
here (■ • ■) denotes the ensemble average. In the thermo- 
dynamic limit, these cumulants saturate to 2/3 in a long- 
range ordered phase, and vanish in a disordered phase. 
For both cumulants, the curves for different L all cross at 
a unique temperature T (see Figs. |2Ia) and (b)), which 
we identify as the transition temperature Tc — 0.65(1) to 
the columnar phase. 

The specific heat per site Cy/N displays a peak, which 
does not diverge in the thermodynamic limit [Fig. [2Ic)]. 
This points towards a second order transition with a crit- 
ical exponent a < 0, or a KT transition. The peak of 
Cy/N is located slightly below the value of Tc determined 
by the cumulants. A featureless specific heat at Tc is a 
strong indication of a KT transition. 

High temperature phase — For T > Tc we calculate the 
dimer-dimer correlation function 

G(x) = (n.(r)n.(r + x)) - 1/16, (3) 

where n.(r) — 1 for a horizontal dimer at site r, and 

otherwise. The constant 1/16 stands for the dimer 
density squared. Taking x = (x, 0) we find that G{x) 
is staggered with x at any T > Tc- We also calculate 
the monomer- monomer correlation function Af(x), pro- 
portional to the number of configurations with two test 
monomers separated by x, and normalized by Af(l) = 

1 Hill [13. M(x) can be computed in the loop build- 
ing process of the MC simulation For a bipartite 
system, M(x) = for monomers on the same sublattice. 

Anticipating the results, we look for algebraic decay of 
these correlators for T > Tc and define exponents and 
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FIG. 3: (Color online) Dimer-Dimer G(x) (left) and 
Monomer-Monomer M{x) (right) correlation functions versus 
distance a; in a log-log scale for various temperatures. 

am by G{x) ~ (— )^x~"'' and AI{x) ^ x^""^ for large x. 
At T = oo, we have the exact results ad = 2 and am = 
1/2 In Fig. we show G{x) and M{x) calculated 
on an L ~ 512 sample for various temperatures T > 
Tc. Both correlators indeed show power-law decay for all 
T > Tc, and the exponents ad and am appear to vary 
continuously with T. The power-law decay is eventually 
truncated around L/2 due to the PBC, making a high 
precision fit for ad{T) and amiT) difficult. 

We therefore turn to TM calculations, which are per- 
formed on semi-infinite cylinders of even circumference 
L. By convention, the space coordinate x is horizontal 
and time t runs upwards. Placing arrows on all lattice 
edges, such that even (odd) sites are sources (sinks) of 
four arrows, we define a (staggered) reference configu- 
ration Co by placing dimers on up-pointing arrows only. 
Superposing a given configuration c with cq gives a tran- 
sition graph with a conserved number Wx of (oriented) 
time-like strings. More precisely, Wx = X^a: ipli^T^) ~ 
nj'{x,t)) e {—"§1 •■•;"§} is independent of t and labels 
the blocks of the TM. The block Wx ^ corresponds 
to a defect of \Wx\ monomers on the same sublattice at 
t — ~oo, and \ Wx\ monomers on the opposite sublattice 
at t = cxj. This can be seen by taking a Wx = configu- 
ration and shifting its dimers along some time-like string: 
this changes Wx by ±1 and introduces a monomer at ei- 
ther end of the string. 

Ordering the TM eigenvalues as AJ^^ > A^"" > . . ., we 
expect A^ (no monomers) to dominate all other eigenval- 
ues. Using standard predictions of conformal field the- 
ory , the L-dependence of A? gives the central charge 
c; the ratio Aj'/A^ determines ad and Xi/X\ gives am- 
Note that the eigenvectors corresponding to (resp. A2) 
are even (resp. odd) under the shift x ^ x + I. 

To account for the interactions in Eq. (Q, the TM 
states must encode vertical and horizontal dimer occu- 
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FIG. 4: (Color online) Dimer-dimer ad and monomer- 
monomer Om, exponents (and their product) obtained by TM 
calculations, and coupling constant g of the CG as obtained 
by MC, all versus temperature T. Inset: central charge c 
versus T. The dashed lines denote the MC estimate of Tc. 



pancies within a row We use sparse-matrix fac- 

torization, and the state space is built from a reference 
state in the given W^j-sector by means of hashing tech- 
niques. Computations were performed up to L = 18, 
with 1 728 292 states in the = sector. 

In Fig.^ we show am and ad as functions of tempera- 
ture (error bars are determined from the finite size scaling 
analysis), increases from 1/2 at T = cxa to 4.0(2) at 
the transition temperature T^, and then diverges below 
Tc (not shown in Fig. ad decreases from 2 at T = cxo 
to 0.26(4) at Tc, and vanishes below (not shown). We also 
plot in Fig. 0] the product am ■ ad, which equals 1 within 
error bars for any T > Tc- The continuous variation of 
both am and ad as functions of T indicates criticality for 
Tc < T < 00. This is corroborated by the central charge 
c as a function of T (see inset of Fig.0J, which is 1.00(1) 
in the critical phase and drops abruptly below Tc. Tc can 
be estimated from the TM by am(Tc) = 4, the jump in c, 
or ac;(Tc) = 1/4. The exact values quoted are accounted 
for by the CG formalism discussed below. The first de- 
termination is the most precise and gives Tc = 0.65(1), 
in full agreement with the MC results. 

Coulomb gas representation — We now account for the 
above findings (KT transition, critical phase with contin- 
uously varying exp onents) by dressing the model's CG 
representation [ll|. Close-packed dimers on the square 
lattice admit a height representation where a height 
z is assigned (up to an overall constant) to each plaque- 
tte in the following way: going counterclockwise around 
a site on the even sublattice, z changes by -f3/4 when 
crossing a dimer and by —1/4 across an empty bond. In 
these units, a monomer corresponds to a dislocation of 



1 in the height. To describe the long- wavelength physics 
of our model, we follow Ref. |Q and coarse-grain the local 
height variable z to a continuum field h{r), for which we 
postulate the following action: 



S ^ j dr [it g\\/ h{r)\^ + V cos {2T:ph{r))] 



(4) 



The first term involving the CG coupling g j2^ accounts 
for the dimer entropy associated with a given coarse- 
grained height /i(r); it dominates in the high-T rough 
phase. The second term is a locking potential of strength 
I^ > (the actual value of V is immaterial) that favors 
dimer ordering in one of the p = 4 "ideal" states (see 
Fig.d satisfying h[v) = 1/8, 3/8, 5/8 and 7/8 (mod 1); 
it dominates in the low-T fiat phase. When T ^ 
the locking potential becomes irrelevant and the height 
model undergoes a roughening transition. 

At T > Tc the height field is thus free and described 
by a Coulomb gas with electric charges e only • Dual, 
magnetic charges m correspond to height dislocations, 
which can be inserted "by hand" through monomers or 
appropriate boundary conditions (cf. the above discus- 
sion of the TM). Units are such that e and m are integers. 
The exponent of an electromagneti c op erator of charge 
{e,m) reads a(e, to) =g^^e^+gm? |llj . We identify the 
leading dimer and monomer correlations by Ud ~ a (1,0) 
and am — a(0, 1). The exact results of Ref. 3 then fix 
f; = 1/2 at T = 00 |^; in general g — g(T) must depend 
on T to account for the continuously varying exponents. 

More generally, electric charges define the vertex oper- 
ators I4(r) = :exp (27rie/i(r)) : appearing in the Fourier 
expansion of any operator which is periodic in h{v). Note 
that Vi appears as a continuum limit contribution to the 
dimer operator 21]. The locking potential corresponds 
to V4 whose marginality determines g{Tc) = 4. 

To check these CG predictions, we calculate g{T) from 
the TM results (see Fig. 0J), using a„i — g. Alterna- 
tively, g{T) is accessible from MC calculations of the 
fluctuations of the winding number {W^) = \{W'^ + Wy). 
The winding number Wx (resp. Wy) is the local height 
difference accumulated by a path winding around the 
torus with PBC in the x (resp. y) direction. (W"^) is 
simply related to g in the free-field point of view |22l |. 
By separating the average tilt from the spatial fluctu- 
ations of the height, the partition function is written 
as a sum over Wx and Wy. The associated Boltzmann 
weight is given by the action of Eq. Q evaluated with 
its classical solution (A/i(r) — 0) given by a linear height 
h{x,y) = -^Wx + -^Wy (and V = 0). Then 



{W^ 



-girn 



(5) 



Note that at g = 1/2, this is a non-trivial statement 
about the Pfaffian formulation of the partition func- 
tion 0. In Fig. 12 we show the MC results for g obtained 
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from {W"^) and Eq. lO, which match with high precision 
the TM evaluation. This validates the CG scenario for 
the transition in our model. 

The CG description also holds for v — +1, i.e., when 
dimer alignment is penalized. In this case, the CG cou- 
pling decreases from 5 — 1/2 at T = cxo to 5 = at 
T — Tf. As T ^ Tt the height model becomes unsta- 
ble and undergoes a tilting transition. When T < Tt 
the winding (W'^) becomes extensive. Analysis of the 
L-dependence of the effective central charge gives the es- 
timate Tt = 0.449(1). More details will appear in Ref.IH 

Doping with monomers — If we allow monomers in 
the model, the height mapping is lost but the CG de- 
scription stays valid. For T > Tc introducing monomers 
is a relevant perturbation [am < 4) and destroys the 
critical phase, but it becomes marginal at Tc- Denot- 
ing the weight of a monomer by e^/-^ and the chemi- 
cal potential by /i, we therefore expect a further criti- 
cal line ^i{T) to emanate from niTc) = —00 and extend 
towards decreasing T. Physically this line separates a 
(weakly) doped dimer crystal from a liquid of dimers 
and monomers. The ground-state has only monomers 
for fj,/\v\ > 1/2 but is a 4- fold degenerate monomer- free 
crystal when /i/|w| < 1/2. The transition is thus first 
order at T = and presumably remains so for small 
T as both phases have gapped elementary excitations at 
/i/|w| = 1/2. However, the transition at low monomer 
doping is second order, so it must cease to be continuous 
at some finite (T* , M* ) • Using a modified TM we estimate 
n = 0.39(4) and ^i^/\v\ = 0.25(10). Along the line, for 
T^, < T < Tc, the TM results are consistent with c = 1 
and a T-independent exponent ad = 1/4. Less relevant 
operators however have T-dependent exponents. 

Discussion and conclusion — In conclusion, we have 
studied a model of dimers on the square lattice with an 
interaction that favors dimer alignment. We find that 
this model undergoes a KT transition separating a high-T 
critical phase with continuously varying exponents from 
a low-T crystalline phase. The transition is understood 
as the roughening transition in the height representa- 
tion or equivalently as a proliferation of e = 4 electric 
charges in the CG framework. Our model captures nat- 
urally all the contents of the CG: electric (dimers) and 
magnetic (monomers) charges, and a T-dependent cou- 
pling constant g driving the transition. Concerning quan- 
tum systems, the rough phase should survive in a large 
part of the finite-temperature phase diagram of the QDM 
(above the melting temperatures of the ordered phases) 
and especially so in the vicinity of the Rokhsar-Kivelson 
point 13 • We thus expect continuously varying exponents 
parametrized by a CG coupling constant at sufficiently 
high T. Our results also indicate that doping the QDM 
would immediately destroy this critical phase. Detailed 
finite-temperature properties of the QDM yet remain to 
be explored. It would be interesting as well to study 
quantum Hamiltonians with ground state corresponding 



to Eq. (UJ 
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